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Abstract. Bose-Einstein condensates are studied in a potential of finite depth 
which supports both bound and quasi-bound states. This potential, which 
is harmonic for small radii and decays as a Gaussian for large radii, models 
experimentally relevant optical traps. The nonlinearity, which is proportional 
to both the number of atoms and the interaction strength, can transform bound 
states into quasi-bound ones. The latter have a finite lifetime due to tunnelling 
through the barriers at the borders of the well. We predict the lifetime and 
stability properties for repulsive and attractive condensates in one, two, and 
three dimensions, for both the ground state and excited soliton and vortex states. 
We show, via a combination of the variational and WKB approximations, that 
macroscopic quantum tunnelling in such systems can be observed on time scales 
of 10 milliseconds to 10 seconds. 
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1. Introduction 

The tunnelling of a particle through a potential barrier is a fundamental effect in 
quantum mechanics 1 . Macroscopic quantum tunnelling is the tunnelling of a many- 
body wavefunction through a potential barrier, and therefore provides a more stringent 
test of the validity of quantum mechanics than the one particle case 0- One place 
where the study of macroscopic quantum tunnelling is particularly experimentally 
accessible is in the tunnelling of a Bose-Einstein condensate (BEC) out of an optical 
trap. Recently, we showed that optical traps, which are of finite depth, can support 
both bound and quasi-bound states, and that the nonlinearity of the BEC mean field 
can be used to tune between them 0QJ|5]. In this article, we calculate the lifetime 
of such quasi-bound states for both the ground state and excited soliton and vortex 
states. This provides a straightforward experimental observable for the occurrence of 
macroscopic quantum tunnelling of a BEC. 

BEC's exhibit different kinds of tunnelling phenomena. We consider the most 
direct generalization of the single-particle case PP, the tunnelling of the mean 
field through a barrier via the Gross-Pitacvskii, or nonlinear Schrodinger equation 
(NLS) [Hj- We emphasize that this is a nonlinear tunnelling problem in the mean- field 
approximation [7| . Experimentally, tunnelling of the mean field in a double- well |5] 
or sinusoidal lattice potential |1U| in configuration space, as well as in spin space for 
multiple-spin-component BEC's ^Ui nas been investigated. There have been various 
theoretical studies related to these experiments [El HS1 HH HS1 HSl 113 CHI EH This is 
to be contrasted with tunnelling of the whole condensate in a variational parameter 
space, as was previously considered in the context of the collapse of a metastable 
attractive BEC in three dimensions [201 EH I—* 1^31 121] and in quantum evaporation 
of a bright soliton in an expulsive harmonic trap |7]. 

In this work, we consider a spherically symmetric trap in the form of a parabolic 
potential times a Gaussian, where the width of the Gaussian envelope is much 
greater than the harmonic oscillator length. This models the optical trap used 
in many experiments on BEC's (see, for instance, [53]). A potential offset Vq at 
the origin, which is important [3| in preventing collapse of attractive condensates 
in three dimensions |26j . can be added with an additional blue- or red-detuned 
laser beam focused at the center of the trap. We take advantage of semi-classical 
methods and employ a variational- WKB formalism to calculate the lifetime of the 
condensate held in such a potential 0. In addition, we calculate three critical values 
of the nonlinearity which could be observed in experiments: the point at which the 
condensate collapses, termed the collapse nonlinearity] the point at which the quasi- 
bound states are transformed into bound states, called the critical nonlinearity, and 
the point at which the condensate is pushed out over the top of the potential and 
quasi-bound states cease to exist at all, called the maximum nonlinearity. The collapse 
nonlinearity is observable as the point where the condensate first implodes. The critical 
nonlinearity is observable as the point where the lifetime first becomes longer than 
inverse experimental loss rates. The maximum nonlinearity is observable as the point 
where the condensate commences to expand. 
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We emphasize that even at the level of the mean field approximation, there are 
features that are distinctly different from the tunnelling of a single particle. As the 
tunnelling rate depends on the number of atoms remaining in the well, the lifetime 
is not simply inversely proportional to the rate as in the linear Schrodingcr equation, 
but must be determined by an integral over each step of the tunnelling. The final 
state of the tunnelling process need not be that of zero probability of there being an 
atom remaining in the well. Instead, a quasi-bound repulsive condensate can decay 
towards a final bound state in which a finite number of atoms remain in the well. A 
quasi-bound attractive condensate can decay towards an unbound state so that the 
atoms spill out over the top of the well. 

In Sec. [3 our application of the variational- WKB method to this nonlinear 
tunnelling problem is explained in detail. In Sec. El the lifetime and the three above- 
mentioned critical points are studied for the ground state in one, two, and three 
dimensions. Excited states are treated in Sec. 14.11 in one dimension for dark solitons 
and bright twisted solitons, and in Sec. l4.2l for vortices in two dimensions. In all cases 
we consider both repulsive and attractive condensates. 



2. Fundamental equations and methods 

The time-independent isotropic NLS with an external potential of the form described 
in Sec. may be written as 

h 2 

-— v 2 * + v^)* + ^|*| 2 * = a*, (i) 

V(r) = (Vo + ±Mu, 2 r 2 ^j exp (-r~ 2 /2£ 2 Gauss ) , (2) 

p OO 

[3 D drr D - 1 \*\ 2 = l, (3) 



where the tildes indicate that the respective variables and parameters are measured 
in physical units. Here ^ is the mean-field wavefunction with the integral norm set to 
unity, M is the atomic mass, ft is a complex eigenvalue which we term the chemical 
potential, ^Gauss is the width of the Gaussian envelope, uj is the angular frequency of 
the parabolic trap, and Vo is the potential offset at the origin. The constant 

/3 li2i3 = 2,27r,47r (4) 

for the spatial dimension D = 1,2,3. The coupling constants are renormalizcd 
appropriately for the dimensionality: 



C/ 1>2 , 3 = (2hu, x )aN , | \l j aN , aN , (5) 

where a is the s-wave scattering length and N is the number of atoms. The transverse 
oscillator frequencies u>± , oj z must be sufficiently high so as to reduce the effective 
dimensionality of the mean field of the BEC [771 EM EM I3U| . 

The NLS can be conveniently rescaled to dimcnsionlcss form: 

- -V 2 * + V(r)* + J/dI^I 2 * = m*, (6) 
V(r) = [Vo + ir 2 ^ exp (-ar 2 ) , (7) 
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Pd drr 15 " 1 ^ 2 = 1, (8) 
Jo 

where all variables arc in units of the harmonic oscillator energy Too and length 
fho = y^h/muj. The coupling constants become 

Ui,2,s = aN/£ ho , aN/£ ho , aN/t ho . (9) 

The parameters Vo and 

a = (4o/^Gau SS ) 2 (10) 

characterize the structure of the potential. For a broad Gaussian envelope, which 
pertains to the experimentally available optical trap described in Sec.^ ct< 1. 

Approximate solutions to Eq. (JSJ) with potential can be obtained via the 
variational method. The ansatz we use in each of the cases to be considered 
below contains a Gaussian factor of form A exp[— r 2 /(2p 2 )], with the amplitude A 
and the width p taken as variational parameters. Substituting the ansatz into the 
normalization condition (jSJ and the Lagrangian of the NLS (for the time being, Im(/x) 
is disregarded), one minimizes the Lagrangian with respect to p and A. Then one 
obtains a system of equations for the nonlincarity Ud and the real part of the chemical 
potential Re(p) in terms of given parameters of the system, D, a, and Vo- The solution 
is stable in the framework of the time-dependent radially symmetric NLS if 

dRc(p)/d\U D \<0, (11) 

which is known as the Vakhitov-Kolokolov (VK) criterion [21] . An important point is 
that our choice of ansatz truncates the solution space by requiring the wavefunction 
decay as r — > oo, eventually enabling us to find quasi-bound states in the form of 
cigenstates with complex eigenvalues. 

To find the imaginary part of the chemical potential one can separately apply the 
WKB approximation ^| E2J • The tunnelling rate 7 in D dimensions is given by the 
standard expressions 

7 =^cxp^-2^ ' dr\ P (r)\^j , (12) 

p(r) = ^2[^-Vf(r)}, (13) 

s 4 T * (14) 
Jo \p(r)\ 

where the endpoints r = 7*1 and r = r 2 are found from setting the semiclassical 
momentum p(r) = and 

Vftf^Vtf+UifaW (15) 



T/ 2 off (r)^(r)--^ + ^|0 2 | 2 , (16) 



T/ 3 off (r)^(r) + ^|0 3 | 2 . (17) 



Here v is the semiclassical oscillation frequency in the well. The imaginary part of the 
chemical potential is then given by 

Im( / i)=- 7 /2. (18) 
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Note that 7 must be multiplied by 2 in the ID case, to account for tunnelling from 
both sides of the well. The transformed wavefunctions in Eqs. (|15fl - i|17|) are given by 
the standard expressions 

0i,2,3 = *, Vr* , r# . (19) 

In the case of axially symmetric excited states in 2D, i.e., vortices, it is useful to 
make the phase winding number explicit: 

02 = 4>2, m e lme . (20) 

Then the transformed wavefunction 02, m leads to the simplified NLS 

1 d 2 , 

where the effective potential is 



V 2.m 02,m = M 02,m , (21) 



^3.m =V X r ) + ^2 + — l02,m| ■ (22) 

Note that the centrifugal barrier changes sign for m ^ 0, so that the semiclassical 
period of oscillation in the well must be redefined as 

dr 

(23) 



W)\ ' 

where tq is the inner turning point created by the centrifugal barrier. The outer 
potential barrier is delimited by 77, r-i. 

The main experimental observable we consider is the lifetime of the trapped 
condensate. This is not given by 1/7, as in a linear system, but must be found 
from the rate equation 

dN/dt = —j(N)N , (24) 

where 7 (AT) is the nonlinear tunnelling rate and N is the number of atoms in the BEC 
remaining in the trap, i.e., between the classical turning points. This leads to the 
integral 

T= / dN/[- 7 (N)N], (25) 

J N 

where T is the lifetime and Nq is the initial number of atoms in the well. Here we have 
assumed that the atomic interaction strength a and the harmonic oscillator length ^ho 
are constant in time. Therefore, in practice, we integrate over the norm Ud, which is 
proportional to N . 

Besides the lifetime, there are three important quantities we will calculate. The 
maximum nonlincarity £/g Lax is the point at which the condensate spills out over the 
top of the well and a quasi-bound state ceases to exist. The critical nonlinearity [7™* 
is the point at which the real part of the chemical potential changes from positive to 
negative, the imaginary part goes to zero, and the quasi-bound state is transformed 
into a genuine bound state. The collapse nonlinearity Uf^ is the point at which the 
condensate collapses in two and three dimensions, as we will explain in more detail in 
the following section. 

There is an important subtlety in the limits of integration in Eq. 1)25(1 . We can 
define an iV crit oc Ug^ and 7V max oc l^ ax . If > 0, then the nonlinearity is always 
repulsive for tunnelling processes. The upper limit of Eq. (|25|l must be replaced with 
(N — N crit )/e + N CI it, and we define the lifetime with respect to decay to a stable state. 
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That is to say, no more than (7V — N clit ) particles will ever leave the well. On the 
other hand, if Up ax < 0, then the nonlincarity is always attractive. The upper limit of 
Eq. must then be replaced with (N — N max ) / e + N max , and we define the lifetime 
with respect to decay to a non-stationary state. In this case, after (jVo — N max )/e 
particles have left the well via quantum mechanical tunnelling, the rest simply flow 
over the top classically. These cases are extremely different from tunnelling of a single 
atom in the linear Schrodinger equation, where the final state always corresponds to 
zero atoms remaining in the well. 

3. Tunnelling of the Ground State in One, Two, and Three Dimensions 

Taking the variational ansatz as a Gaussian, 



one minimizes the Lagrangian corresponding to Eq. ©. The result is a system of 
Euler-Lagrange equations for the chemical potential and nonlinearity [221 El ■ These 
take the form 



where T(D/2) is a Gamma function and D is the dimensionality. 

The transition from a quasi-bound to a bound state occurs when the real part 
of the chemical potential changes sign from positive to negative, since our potential 
V(r) — > + as r — > oo. The nonlinearity at which this occurs we termed £/™'. In Fig.Q] 
is shown the dependence of C/™' on the potential offset Vb, as defined in Eq. JJJ. The 
cases of I? = 1,2,3 dimensions are illustrated in separate panels in the figure. The 
different curves in each panel pertain to different values of the trap-shape parameter 
a. These two continuous parameters, Vo and a, and the discrete one, D, are the only 
given constants in the system. The potential offset Vq can be controlled by a red- or 
blue-detuned laser focused at the center of the trap, as previously mentioned. 

However, the existence of the transition from a quasi-bound to a bound state does 
not mean that the ground state solution we obtained is stable. When the nonlincarity 
is negative, so that the atoms attract each other, the condensate can collapse. In three 
dimensions with V(r) = 0, any initial condition leads to collapse (201 . As is well-known 
for BEC's, the imposition of a nonzero potential V(r) can create a metastable solution 
with a lifetime much longer than that of the BEC, so that it is experimentally stable. 
This requires a nonlinearity which is not too strongly negative, U3 > U^° n , where C/3 011 
is the critical point for collapse. This critical point can be determined variationally 
or numerically |35l 1361 12] • In order to observe the transition from a quasi-bound to a 
bound state, it is necessary that £/|° n < U§ rlt , as we showed recently [2]- The critical 
nonlinearity for collapse corresponds to an eigenvalue pf° n , while [7™* corresponds to 
(i = 0. Therefore a simple way to state this condition is 




(26) 



Re( M ) = \ {(-4 + D)p-^ + (1 + ap 2 )- 2 -^ 2 



x [4V + (4 + D)p 2 ~a(D + 4Vb«)p 4 ] } , (27) 

U D = Y(D/2)2- 2+D l 2 p p - 2+D {\ + ap 2 )- 2 - ' 2 

x {-2(1 + ap 2 ) 2+D ' 2 + p 4 [2 - AV a - a(D + W a)p 2 }} ,(28) 



coll 



< 0. 



(29) 
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In Fig.|2», ^ co11 is shown as a function of the trap parameter a for Vq — 0. Clearly, 
in order to observe the transition from a quasi-bound to a bound state in three 
dimensions, one needs a finite offset Vq < 0, as was shown in Ref. [3]. For a<Cl, the 
effect of a nonzero offset Vb on p c ° n is essentially linear: p c ° n ~ /^ ooll (Vb = 0) + Vq. 
For instance, to observe the transition from the quasi-bound state to the bound one, 
it is required that Vb < —0.273 for a = 1/4. For larger a, the relation between p, and 
Vq becomes nonlinear. For example, in the case of a = 1, the necessary condition is 
Vb < —1.495. Since such values of a are experimentally irrelevant, we do not illustrate 
this regime. 

Inspection of Fig. Ufa) reveals that the curve terminates at a = 0.2. For a > 0.2 
and Vo = 0, there is no quasi-bound state at all. Another way of stating this is that 
for any a, there is a critical Vq "* such that 

V < V CTit (30) 

is required to obtain any stationary solution of the form described by the ansatz 
of Eq. (|26|) . A plot expressing this condition is displayed in Fig. Gib). This is an 
important consideration in choosing the correct experimental regime to observe quasi- 
bound states in three dimensions. 

In ID there is no collapse at all. In 2D and in free space, initial conditions lead to 
collapse for U% < J7|° n and expansion for U2 > U^ 011 , while for U2 = i7f oU an unstable 
stationary state is obtained, known as the Townes soliton This weak collapse is 
different from the strong collapse which occurs in 3D |2fij . An appropriate external 
potential can stabilize the expanding regime, producing a stable solution, rather than 
a metastable one as in 3D (see, for instance, the discussion in Ref. 001 )• The value 
[/2 011 which determines the collapse threshold can be derived from Eqs. (|27[) - i|28[) by 
taking p — > + . Then p — > —00 and 

{7| oll = -27r. (31) 

A third critical nonlinearity is given by the point at which the eigenvalue p = 
max[V A (r)], i.e., where the quasi-bound state ceases to exist. We called this maximum 
nonlinearity £/™ ax . The range of allowed nonlinearities for obtaining a stationary state 
is therefore given by {J™ 11 < U2 < C/" lax - I n this range all solutions are stable or quasi- 
stable. Outside of this range there are no stationary states. Figure illustrates the 
dependence of C/™ ax on the other parameters in the problem. The collapse point is 
shown as a dashed blue line. Note that the inequalities [Z™ 11 < C/| rlt < U™ ax hold for 
all shapes of the trap, unlike the 3D case. 

A technical point concerning two dimensions is that there is always an additional 
solution with large p which is radially unstable within the context of our Gaussian 
ansatz and the VK criterion. This corresponds to a ring of condensate around the 
trap. Its variational form and stability can be more appropriately investigated with a 
vortex-like ansatz, as we briefly discuss in Sec. 14.21 

A clear experimental signature for a quasi-bound to bound state transition is a 
change in the lifetime of the condensate. Suppose that the lifetime due to tunnelling 
of the wave function is on the order of or much less than that imposed by three body 
processes, scattering with the background gas in the vacuum, and other sources of 
loss. One can then determine the difference between a quasi-bound and bound state 
by measuring the number of atoms in the trap as a function of time. It is necessary 
to choose the right experimental parameters in order to have a realistically observable 
lifetime. For instance, typical BEC lifetimes are on the order of ten to a hundred 
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seconds. However, the mean- field induced dynamical time scales are typically on 
the order of milliseconds. Therefore, a lifetime due to tunnelling on the order of 10 
milliseconds to 1 second is desirable. 

In Fig. we show the lifetime for several values of a. Panels (a)-(c) pertain to 
one, two, and three dimensions. The lifetime is calculated from Eq. as discussed 
in Sec. El The leftmost endpoint of the curves corresponds to t/fj 1 *. When £/|> rlt > 
the final state is bound and has a nonzero number of atoms. Three examples are 
shown for D = 3 in Fig. 0{c) . The rightmost endpoint of the curves corresponds to 
(/max \Yhen C/™ x < 0, the final state is unbound and the remaining atoms spill out 
over the top of the trap without tunnelling. Two examples are shown for D = 2 in 
Fig. lib). 

Since the lifetime is scaled to the trap frequency uj 1 one can convert the range of 
T to milliseconds easily. For instance, for a trap of angular frequency u) = 2ir x 100 
Hz, the range T = 10 to 10 4 shown in Fig. 21 corresponds to 16 ms to 16 s. Thus the 
trap parameters chosen for the figure result in experimentally observable lifetimes. 

The relationship between the tunnelling rate, defined by Eq. (|12f) . and the lifetime, 
defined by Eq. Ij25(l . can be very non- intuitive with respect to what is known from 
single-particle tunnelling in the linear Schrodingcr equation. For example, in Fig. 01c) 
the lifetime approaches zero on the left hand side as U3 — > f/| rlt > 0. This can 
be physically understood in terms of the tunnelling of the wavefunction through the 
effective potential. As the real part of the chemical potential approaches zero from 
above, so that the quasi-bound state becomes a bound state, the attractive nonlincarity 
pulls down the potential barrier. That is, the effective potential barrier shrinks. In 
the case of three dimensions, and for Re(/x) < 0, it is well-known that the attractive 
nonlincarity can dominate over both the kinetic energy and the potential energy, 
leading to strong self-focusing and collapse. Here, one observes that the attractive 
nonlincarity already dominates over the potential energy barrier for Rc(^) > 0. This 
causes the effective potential barrier to disappear for a very narrow region of U3 near 
the formation of the bound state, so that the lifetime goes to zero. This region is 
so narrow that it is not experimentally relevant, since number fluctuations cannot be 
finely controlled in BEC experiments, and Us oc N. 

One can more rigorously understand the limiting values of the lifetime in the form 
of the power law of the tunnelling rate 7 as U3 — > (t/| rlt ) + ■ Suppose the power law 
takes the form 

j(U D ) ~ (U D - Ul^Y , (32) 

where p is a real positive constant and [/|, nal can be zero, C/^ ax , or C/£) rlt , depending 
on the allowed domain of quasi-bound states (see the discussion of the appropriate 
integration limits following Eq. 1)25(1 ^1. By expanding the integral according to 
e dxf(x) ~ e/(a), one finds that for < p < 1 the lifetime approaches zero, 
despite the rate approaching zero. For p > 1 the lifetime approaches infinity, while 
for p = 1 it approaches a nonzero constant. We find numerically for the ground state 
that whenever C/£ nal = Ug^ > or U% naX = Uff^ < the power law of the rate is 
such that the lifetime approaches zero at these critical points. For C/f, nal = 0, which 
occurs when < and Up ax > 0, p — 1 and the lifetime approaches a nonzero 

constant. 

Finally, we note one other important point about the lifetime. For Ud = 0, where 
there is no mean field, one recovers the correct single-particle T = 1/7, as can also be 
seen by consideration of Eq. Ij25(l . This is a limiting case in the center of a smooth 
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lifetime curve in Fig.0](see also Figs. HJ] and [5] below). Therefore, it strongly supports 
the validity of the mean field approximation. 

4. Tunnelling of Excited States 

4-1- Soliton Tunnelling in One Dimension 
The variational ansatz 

<& sol = Arexp(-^\ (33) 

models the excited state with a single node in one dimension. This is a dark soliton 
for U\ > 0, i.e., a repulsive BEC in ID, and an antisymmetric, or bright twisted soliton 
for U\ < 0, i.e., an attractive BEC in ID. The variational equations for the chemical 
potential and nonlincarity may be derived by means of the methods outlined in Sec. [3 

Rc(m) = \ {-9p- 2 + (1 + ap 2 )- 7 / 2 [AVo + (15 - l6V a)p 2 

-5a(3 + 4V *)p*] } , (34) 

(-2 + (1 + <V)- 7 /V [2 - W a 
P L 

-a(3 + W Q a)p 2 ] } , (35) 

A plot of (7f rlt , the dependence of the critical nonlincarity for the transition from a 
quasi-bound state to a bound state, is shown in Fig. [5j The transition point is only 
very weakly dependent on a for a«l. This regime is an experimental relevant one. 

Figure [fJl shows the main experimental observable characteristic of quasi-bound 
solitons, namely, the lifetime. The node causes the wavefunction to become somewhat 
broader than the ground state. This requires smaller values of a in order to provide 
sufficiently thick walls to hold the quasi-bound condensate, in comparison to Fig.^a). 
As was true for the ground state, the leftmost endpoint of the curves corresponds to 
U" A , while the rightmost endpoint corresponds to [7{ nax . For a = 1/6 and Vq = 0, 
one finds f7™ ax < 0. Thus, after losing some atoms via tunnelling through the barrier, 
the final state takes the form of two bright solitons with a phase difference of it which 
travel in opposite directions away from the potential center. 

4-2. Vortex Tunnelling in Two Dimensions 

Under the assumption of an axisymmetric stationary state in two dimensions, one may 
take the variational ansatz for a vortex state as 

2p 2 

where m is the winding number, or topological charge. Then, transforming the ansatz 
according to Eq. (|20|l and minimizing the Lagrangian with respect to Eq. I|21|) with 
effective potential (|22|l , one derives the following system of equations for the chemical 
potential: 

r*oo - ^ttV^ \- r -^ + (i + «P 2 r 3 -™ 

2T(1 + m) I p 2 - 

x [-2V (1 + ap 2 ){-\ + (1 + 2m)ap 2 )r(l + m) 



A r m cxp ( -— ) exp(im(9) , (36) 
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-p 2 (-3 + (1 + 2m)ap 2 )r(2 + m)]} , (37) 

U2 = _ r( 1 + ^r (2++2: )^ {l + (l + apV _„ 

xp 4 [-1 + 2F a + a(l + m + 2U a)p 2 ] } , (38) 
As a consistency check, note that, for the case m = 0, one recovers the same result as 
given by Eqs. and <EEJ> wi *h O = 2. 

The collapse point may be obtained by expansion of Eqs. I|37|l and (|38fl in small 
p, in the same way as was done in Sec. |3J One thus finds 

u _ 2^(2 + m) 
U * -~ r(l/2 + m) ■ (39) 

For m = 0,1,2,3 one obtains J7 2 coU = -2tt, -8tt, -1&k, -(128/5)tt. The presence of 
the vortex decreases the collapse point, as was pointed out some time ago [351 1371 13"%] . 
The critical point for a transition from a quasi-bound to a bound state is shown in 
Fig. for m — 1,2,3. Like the ground state in two dimensions, the nonlinearity is 
constrained to the range C/™ 11 < U2 < U™ ax . In Fig. [5] the allowed range is shown as 
a function of a and Vq for m = 1,2. 

The lifetime of the condensate with a vortex is shown in Fig. [5] for m = 1,2. For 
higher winding numbers the ansatz becomes broader, and it is necessary to choose 
smaller a in order to contain the vortex. This is reflected in the range of a in Fig. [51 
a = 1/8 to 1/12 for m = 1 and a = 1/14 to 1/18 for m = 2. Note that, unlike for 
the cases illustrated in Figs. E] and [S] here T — > 00 as U2 — > C/| rlt . In Fig. E{a), for 
a — 1/8 the atoms first tunnel through the barrier. When U2 reaches J7™ ax < 0, the 
condensate spills out over the top of the potential barrier. Then it expands towards 
r = 00 as a bright ring soliton j^Hl EH] , pushed outwards by the Gaussian tail of the 
potential. At a critical value of the ring radius, it undergoes azimuthal modulational 
instability |3T] . 

We have not considered azimuthal instabilities. These should be manifest 
for nonzero winding number in the attractive case in 2D in certain parameter 
regimes [SHIES]- A simple stability criterion is that the wavelength of modulational 
instability be longer than the ring circumference 2irp (for details, see 001 and references 
therein). A more sophisticated treatment via a full linear stability analysis (using the 
Bogoliubov equations) is left for a future work. 

We note that for nonzero winding number there is an additional stationary state 
which has a radius greater than that of the potential peak. For repulsive nonlinearity, 
this is a radially unstable pinned vortex with an initial core size artificially controlled 
via the trapping potential. For attractive nonlinearity, this state can be radially 
stabilized for sufficiently strong nonlinearity, as may be shown via the VK criterion of 
Eq. ([11[) . However, in this regime we expect the azimuthal instability to dominate, so 
we have not presented it here in any detail. 



5. Conclusion 



In this study we suggested a straightforward macroscopic quantum tunnelling 
experiment which requires little or no modification of existing experimental apparatus. 
We assumed a trap consisting of a parabolic potential times a Gaussian envelope, 
which models typical optical traps used in experiments. Such a trap supports both 
bound and quasi-bound states. Using a variational- WKB mean-field formalism, we 
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calculated four experimental observables: the lifetime T of quasi-bound condensates; 
the maximum nonlincarity Up a * for which a quasi-bound state exists; the critical 
nonlinearity [7™* for which a quasi-bound state becomes a bound state; and, for two 
and three dimensions, the collapse nonlinearity U^ 11 < 0. 

By adjusting the initial number of atoms and/or the atomic interaction strength, 
as may be achieved via a Fcshbach resonance [?, 1441 145| . we showed how a quasi- 
bound condensate can be adiabatically transformed into a bound one, and vice versa. 
We presented in detail the parameter regimes in which such a transformation can 
be performed for both repulsive and attractive condensates. In two dimensions, we 
found that the relations ?7™ ax > £/f rlt > C/™ 11 always hold; in three dimensions, it is 
necessary to add a negative offset to the potential at r = in order to observe the 
quasi-bound to bound transition without provoking collapse of the condensate. This 
offset can be constructed with an off-resonant blue- or red-detuned laser focused at 
the origin of the trap. 

We found that, for reasonable experimental parameters, one could observe the 
tunnelling of a BEC through the walls of an optical trap on time scales of 10 
milliseconds to 10 seconds. We showed that tunnelling can lead to a final state 
which is quite different from that obtained via the linear Schrodinger equation. For 
Up > 0, an initially quasi-bound repulsive condensate approaches a bound state 
with a finite number of atoms remaining in the well. For U^ ax < an initially quasi- 
bound attractive condensate approaches an unbound state for which the atoms spill 
out over the top of the well. This is observable as a pair of counter-propagating bright 
solitons in one dimension, and as an azimuthally unstable bright ring soliton in two 
dimensions |3§1E!0"| . 

We showed previously that our variational- WKB method accurately determines 
the tunnelling rate and critical points at the level of 1% or better, except near the 
collapse points, where it is accurate to within about 10% 0. It is also possible to 
test our method via exact solution with a small number of atoms. For example, it 
has been shown that mean field effects are strongly evident in as few as three bosons 
confined in an external potential |46| . 

Finally, we note that there is a direct analogy between the operation of a laser 
and the tunnelling of a BEC through the walls of an optical trap. In a laser operating 
well above threshold, the coherent state in the cavity is robust against removal of 
single photons via tunnelling through a thin barrier. This is the essence of a coherent 
state. The mean field of the Bose condensate in an optical trap is phase-coherent and 
robust in the same way. Therefore we expect our mean field approximation to be an 
excellent one for a large initial number of atoms in the trap. When the initial number 
of atoms becomes very small, the condensate is liable to phase fluctuations. However, 
as our lifetime definition and figures have the right single-particle limit for zero mean 
field, we expect that the mean field theory also gives good results in this region. 
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Figure 1. Ground state transition. Shown is the critical nonlinearity for a 
transition from a bound to a quasi-bound state, as a function of the potential 
offset Vo- (a) One dimension, a = 1,1/2,1/10 (dotted, dashed, solid curves), 
where a = ^jL/^q USB - There is no collapse in one dimension. (b) Two 
dimensions: a = 1,1/2,1/10 (dotted, dot-dashed, solid curves); the critical 
nonlinearity for collapse, f/^ " = —2ir is indicated by the blue dashed curve. 
Note that the condition t/| rlt > f/| o11 always holds, (c) Three dimensions: 
a = 1,1/2,1/3,1/4,1/5,1/6,1/7 from bottom left to top right. The collapse 
point is described in Fig. l2l 
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Figure 2. Ground state parameters, three dimensions. (a) Shown is the 
dependence of the collapse point on the trap-shape parameter a for zero offset 
(V = 0). For nonzero V and a < 1, (U co11 ~ /i coll (Vb =0) + V . In order to obtain 
an observable transition from a quasi-bound to a bound state, it is necessary that 
^coll ^ q ^) Shown is the critical offset V crlt needed to obtain any quasi-bound 
state for a given a: Vq < V max . 
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Figure 3. Ground state parameters, two dimensions. (a) Shown is the 
dependence of the maximum nonlincarity (7™ ax to obtain a quasi-bound state 
on the potential offset Vb- The curves show a = 1/2, 1/3, . . . , 1/10 from bottom 
left to top right. The dashed blue line shows the collapse point, U^ 011 = —2n. In 
this harmonic times Gaussian potential, one is therefore limited to the range of 
nonlinearities £/| o11 < U2 < C/™ ax . Note that for C/ 2 < i/™ 11 there is no stationary 
state at all. 
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Figure 4. Ground state lifetime. Shown is the lifetime T of the BEC, scaled 
to the trap frequency w, for zero potential offset Vq and for several values of the 
trap-shape parameter a = ^ /^Q allss i as a function of the initial nonlinearity 
U D oc N. (a) D = 1, Vb = 0, a = 1/4,1/5,1/6,1/7,1/8 (from bottom to 
top), (b) D = 2, V = 0, a = 1/4,1/5,1/6,1/7,1/8 (from bottom to top), (c) 
D = 3, Vo = -2, a = 1/4,1/6,1/8 (from bottom to top). The leftmost point 
of each curve approaches the value C/^ 1 *, where fi — > + , the quasi-bound state 
disappears, and the lifetime T — > oo. The rightmost point corresponds to C/]^ ax , 
where a quasi-bound state is no longer supported by the potential. 
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Figure 5. Excited state transition, soliton. The critical nonlincarity for a 
transition from a quasibound to a bound state is shown, in analogy to Fig. Eta). 
Positive values of U^ Tlt pertain to the dark soliton; negative values are for a bright 
twisted soliton. The dashed line divides the two regimes. The curves correspond 
to increasing values of a = ^o/^Gauss : a = 1> V^i 1/9, . . . from bottom to top of 
the y-a,yAs intercept. For sufficiently small a, all curves lie on the final dark line. 




Figure 6. Excited state lifetime, soliton. Shown is the lifetime of the excited state 
with a single node in one dimension, i.e., a dark soliton for U\ > and a bright 
twisted soliton for Ui < 0. The parameters are Vo = and ct = 1/6,1/8, 1/10 
from left to right. 
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Figure 7. (color online) Excited state transition, vortex. The critical nonlinearity 
for a transition from a quasibound to a bound state is shown, in analogy to 
Fig-EIb), for values of the winding number m = 1 (solid black curves), m = 2 (dot- 
dashed blue curves), and m = 3 (dot-dot-dashed red curves). Each set of curves 
rises along the y-intcrcept for decreasing a as a = 3~ 2 , 4" 2 , . . . , 1CT 2 . The 
dashed horizontal lines show the nonlinearities for which the solutions collapse: 
[/coll _ _8tt, -16tt, -(128/5)vr for m = 1,2,3. Note that Uf A > C/| o11 for all 
Vo,ct, m. 
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Figure 8. Excited state parameters, vortex. Shown is the dependence of the 
maximum nonlinearity for confinement of a vortex, C/™ ax , on the parameters 
Vo,ct, m. Here Vb is the potential offset, a the trap-shape parameter, and m the 
vortex winding number. Panels (a) and (b) illustrate m = 1 and 2, respectively, 
while a = 1/2, 1/3, 1/10 from bottom to top in each panel. The asymptotic 
right hand value of each curve represents the critical parameter set for which a 
vortex is pushed outside the potential, i.e., where it becomes pinned rather than 
confined. 
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Figure 9. Excited state lifetime, vortex. Shown is the lifetime of the condensate 
with a vortex in two dimensions, (a) Winding number m = 1 and trap-shape 
parameter a = 1/8, 1/10, 1/12. (b) m = 2 and a = 1/14, /, 1/16, 1/18. The 
leftmost point of each curve, which is not shown on the plot, approaches 
where fi — > + , the quasi-bound state disappears, and the lifetime T — > oo. The 
rightmost point corresponds to U™ 11 *, where a quasi-bound state is no longer 
supported by the potential (see Fig. Isl. 



